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Introduction 

Throughout, let A := Z ( "bi- infinite" case) or A := Z>o = {0, 1, . . . } ("semi- 
infinite" case). Consider the set of equations 

^ = A"moo, ^ = -m^(Ay)n^ n = 1, 2, . . . , (0.1) 

on bi- or semi-infinite (i.e., Ax A) matrices moo = i^ooit, s), where the matrix 
A = {Sij^i)ij^A is the shift matrix, and A""" its transpose. In [0, |^, it was 
shown that Borel decomposing]] 

moo(t, s) = {l^ij),j^A = S^'S2, for "generic" t, s G C°°, (0.2) 

into lower- and upper-triangular matrices Si{t, s) and 5*2 (t, s), leads to a two- 
Toda (two-dimensional Toda) system for Li := SiASi'^ and L2 = S2A^ 82^, 
with r-functions given0 by 

Tn(t, s) = detninit, s), n E A. (0.3) 

This paper deals with skew- symmetric initial data moo{0, 0). As readily seen 
from formula (0.1), the 2- Toda fiow then maintains the relation moo{t,s) = 
— moo(— s, —t)~^, and hence, by formula ( p.3| ) and the interpretation of its 
right hand side in footnote ^ 

r„(t,s) = (-l)V„(-s,-t). (0.4) 

The main point of this paper is to study the reduction s = —t, as used in 



the theory of random matrices in H. Peng's doctoral dissertation |]T6| . When 
s —t, formula ( p.4|) shows that in the limit the odd r-functions vanish, 
whereas the even r-functions are determinants of skew-symmetric matrices. 



In particular, the factorization (p.2|) fails; in fact in the limit s —>■ —t, the 



^ Here "t, s G C°°" is an informal way of saying that t and s are two sequences of 
independent scalar variables; a function of those variables may be defined only in an open 
subset of C°° X C°°, or may even be a formal power series in t and s. 

^ This formula will be used mainly in the semi-finite case, with m„ = {iJ,ij)o<i.j<7i- In 
the bi-infinite case, m„ — {^ij)-oo<i.j<n and the determinant is interpreted as 



lim det(^ij)_2fe<i.j<9, (*) 

k — >CXD 



assuming the limit makes sense. 
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system leaves the main stratum to penetrate a deeper stratum in the Borel 
decomposition. In this paper we show this stratum leads to its own system, 
whereas in a forthcoming paper with Horozov we show this system is 
integrable by producing its Lax pair. 

Thus, we are led to considering Pfajjians: 

frXt) := Pfaff m„(t, -t) = (detm„(t, -t)Y/^ = r^it, -tf/\ (0.5) 

for every even n & A (the same remark as in footnote ^ applies here). The 
"Pfaffian f-function" is itself not a 2-Toda r-function, but it ties up remark- 
ably with the 2-Toda r-function r as follow^: 

T,^{t, -t - [a] + m = f2n{t)f2n{t + [«] - [P]) 

r2n+i{t, -t - [a] + [f3]) = {(3- a)f2n{t - [/3])r2„+2(t + N). 

When /3 — > a, we approach the deeper stratum in the Borel decomposition of 
moo in a very specific way. It also shows that the odd r-functions r2n+i(t, —t— 
[a] + [P]) approach zero linearly as /5 — a, at the rate depending on a: 

lim T2n+iit, -t - [a] + [0\)/{(3 - a) = f2„(t - [a])f2n+2{t + [«])• 

Equations ( p.6| ) are crucial in establishing bilinear relations^ for Pfaffian f- 
functions, where n, m G A: 

^ Pj{~2y)e^ -^'^'Pk{-D)f2n ■ r2m+2 

j,k>0 
j-k=-2n+2m+l 

+ Pji2y)e>^-y''''pk{D)f2n+2-f2m = 0. (0.7) 

k-j=-2n+2m-l 

This is a generating function for the Hirota equations satisfied by f{t): for 
each n and m, after expanding ( p.7| ) into a power series in y = {yi,y2, . . .), 
the coefficient of each monomial in y gives a Hirota equation. For example, 

Ma]:=(a,^,4,...) 

^ d = ( J^, ijj, . . . ), and £» = (Di, iDj, ^D^, . . . ) is the corresponding Hi- 
rota symbol, i.e., P(i))/ • g :== P{d/dyi, \d/dy2, ^d/dya, . . . )f{t + y)g{t - y)\y=Q for any 
polynomial P\ and pk are the elementary Schur functions: 'YI'q' Pki^^z^ '■— exp(^j" tiz"^)- 
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for m = n — 1, the coefficient of the hnear terms (2/^+3, to be more specific) 
gives 

Pk+ii-D) - f2n ■ f2n = Pk{D)f2n+2 ' r2n-2, (0.8) 

where k G Z>o and 2n — 2 E A. For A; = 0, this equation can be viewed 
as an inductive expression of r2„+2 in terms of f2n-2 and derivatives of f2n- 
These equations aheady appear in the work of Kac and van de Leur , in 
the context of the DKP hierarchy. On the exact connection, see forthcoming 
work by J. van de Leur UlSj . 

In analogy with the 2-Toda or KP theory, one establishes Fay identities 
for the Pfaff f-functions. In this instance, they involve Pfaffians rather than 
determinants: 



pfg^g f i^j - Zi)f2n-2{t - [Zj] - [Zj]) 
V r2n{t) 



l<i,j<2k 

T2n-2kit - YlTi^i 

= A{z) ^. (0.9) 

T2n{t) 



In the semi-infinite case, the latter has a useful interpretation in terms of 
Pfaffians of "Christoffel-Darboux" kernels of the form 



n-l 



Knifi, A) = e^"*'('''+^') ^(^g2fc(t, A)g2fc+i(t, /u) - q2k{t, fi)q2k+i{t, A)), 

(0.10) 

where the qm{t, A) form a system of skew-orthogonal polynomials [0. This is 
the analogue of the Christoffel-Darboux kernel for orthogonal polynomials. 
So, formula (0.9) can be rewritten as 



PfaS{Kn{zi,Zj))i<ij<2k = Y\ X(t;2;i)r| 



5 



(0.11) 



ordered 



where X(t; z) is a vertex operator for the corresponding Pfaff lattice (see 
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This vertex operator also has the remarkable property that for a Pfafiian 
f-f unction, 

f + aX(A)X(/x)f 

is again a Pfaffian f-function. 

As was shown in ^ , the 2-Toda lattice has four distinct vertex opera- 
tors. Upon using the reduction s = —t, the 2-Toda vertex operators reduce 
to vertex operators for the Pfaff lattice. This enables us to give the action 
of Virasoro generators on Pfaff f-functions, in terms of the restriction (to 
s = —t) of actions on 2-Toda r-functions: 

(jf (t) + (-l)Vf T2n{t,s)l=.t = 2f2n{t)jt\t)r2n{t). 

Finally, we discuss two examples, a first sketchy one, involving a semi-infinite 
Pfaff lattice and matrix integrals; this is extensively discussed in A 
second example, genuinely bi-infinite, will be given in the context of curves 
with fixed point free involutions l, equipped with a line bundle jC having a 
suitable antisymmetry condition with respect to l. 



1 Borel decomposition and the 2-Toda lattice 

In 1^, we considered the following differential equations for the bi-infinite 
or semi-infinite moment matrix rrioo 



A"moo, ^ = -moo(A')^ n=l,2,..., (1.1) 

where the matrix A = (5jj_i)jjg^ is the shift matrix; then (|1.1|) has the 
following solution 

m^(t,s) = e^*"^"moo(0,0)e~^^"(^^)" (1.2) 

in terms of the initial data moo(0, 0). 

Assume moo allows, for "generic" (t, s), the Borel decomposition rrioo = 
S'f for 

g Q _ flower-triangular matrices! 
^ ~ ' (^with I's on the diagonal y 

g Q _ f upper-triangular matrices 1 
^ ~^ ' (^with non-zero diagonal entries J ' 
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with corresponding Lie algebras g^, g^. Assume moreover that moo, 5*1 and 
5*2 are nice in the sense of Remark at the end of this section. Then setting 
Li := S^KS^\ 

g dm^^-i ^ I S,id/dh)iS^'S2)S,' = -S^S^ + S^S,' E g^ + g+, 
' ' [S.A^-m^S,' = S.A^S^' = = + e g. + g+; 

the uniqueness of the decomposition g^ + g^ leads to 
dSi 1 dS2 1 



-s^' = ra-, ^s^' = {LIU. 



Similarly, setting L2 = 82^^ 82'^, we find 

9Si I , , 882 I 



dsr, ds 



This leads to the 2-Toda equations ^Tj for Si, S2 and Li, L2: 



d d 

gr-S{i} = g-S{i^ = ±{L-)^S^iy (1.3) 



^ = [ra+,^.], ^^ = m)^,ui z = i,2,..., (1.4) 

and conversely, reading this argument backwards, we observe that the 2-Toda 
equations ( |1.3| ) imply the time evolutions ( |1.1| ) for moo- 

The pairs of wave and adjoint wave functions ^ = (\E'i,\E'2) and \E'* = 
(^t,^^), defined by 

^|i|(t,s,z) = e^i Uf' 5{i}X(^), 



1 5 



where x{^) is the column vector (z")„g^, satisfy 

= {z, 2"^)^, L*^* = (z, ;z-^)^*. 
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and 



dt, 
_d 
ds, 
_d_ 

dtn 

d 

ds„ 



[1.6) 



which are equivalent to ( |1.3| ), and are further equivalent to the following 
bilinear identities,^ for all m,n E A and t, s, t', s' G C°°: 



/ ^,n{t,S,z)^l„Xt\s: ,Z')^ = (f ^2n{t,S,z)%^{t',S,' ,Z . . 

Jz=oc ^TTZ^; J^^Q 2mz 



[1.7) 



By 2-Toda theory |]T^ the problem is solved in terms of a sequence of 
tau-functions 

Tn{t,s) = detmn{t,s), (1.8) 
with mn(t,s) defined in (and detm„ interpreted as in) footnote |: 

1 (uii(t,s)) ^. (bi-infinite case), 

m„(t,s) := <^ y ^ -~<*'^<" ^ ^ (1.9) 

{Hij{t, s))^^.^.^^, (semi-mfinite case, with tq = 1), 



as 



^,it,s;z)= M-}f}'K ^r''^\- 

V rn{t,s) 

^,(t,s;z) = ( w(M-M) ^Er^.^-'^n 

V Tn{t,s) 

nit,s,z) = ( i!^±i(i±fe:Me-ErM', 

V Tn+l{t,s) 

^*(t s z)- ( !!^i^li±Me-Sr«»^-%-" 



5 



[1.10) 



z ^— n 

I 

neA 



^ The contour integral around z = cx3 is taken clockwise about a small circle around 
z — oo, while the one around z = is taken counter-clockwise about z = 0. 
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Note O and (|nO| ) yield 



h{t, s) := (diagonal part of 5*2) = diagi ^ ) . (l-H) 



Formulas ( |1.7D and ( |1.1(J| ) imply the following bilinear identities 



r, 



„(t - [z-'],s)T^+,it' + [z~'], s')e^^'~'^-'>'z--'^-'dz 



z=oo 



s - [^])r„(t', + [z])e^^^^^-'''^^-'z--^-'dz, (1.12) 

z=0 

where m, n & A, satisfied by and characterizing the 2-Toda r-functions. 
Remark: 



In the bi-infinite case the factorization moo = 81^82 in (p^ ) or the de- 



terminant formula ( p. 3D may fail to make sense. Nevertheless, we can take 
(|0.4| ) as a starting point, use (|0|) to define f up to the sign, and make sense 
of the r-side of the whole story. 

In the bi-infinite case, factorization as in (|0.2|) is not unique in general. 
This is responsible for the Backlund transform of a finite band matrix having 
a continuous family of solutions. However, it means the matrix multiplication 
may not be associative in the bi-infinite case: if 

Q— 1 c c' 

with 5*1, S[ G G- := {lower triangular matrices with I's on the diagonal} 
and 5*2, 52 G := {upper triangular matrices with non-zero diagonal}, and 
if the matrix multiplication was always associative, we should have 

S[S^^ = S'^S^^ eG+nG- = {1}, 

so that 5*1 = S[ and 52 = 6*2. 

The associativity is important in establishing the relation between equa- 
tion ( p.l| ) and the 2-Toda flows on {Si, 82)- Moreover, we are mainly inter- 
ested in the semi-infinite case, in which the associativity clearly holds. So 
we assume that, for generic (t, s), Si, S2 and moo actually belong to suitable 
subgroups G'j^ of G± and a suitable subspace Mat' of the space Mat of all 
infinite matrices, respectively, in which the multiplications 

G'_xMat' ^ Mat' Mat'xG"+ ^ Mat' 

[o,m) 1-^ bm [fn,J) ^ 
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are associative: 

{SS')m = S{S'm), m{S"S"') = {mS")S"' and {Sm)S" = S{mS") 

hold for any S, S' G G'_, S", S'" G G; and m e Mat'. 

For instance these conditions are satisfied if the entries of every 

matrix in G'^ and Mat' tend to quickly enough as i — > —oo uniformly in 
j > z + a, and as j —>■ — oo uniformly in i > j + a, for some constant a. 

2 Two-Toda r-functions versus Pfaffian f-functions 

In this section, we assume either the matrix moo is semi-infinite, or det m„ can 
be interpreted as in formula (*) in footnote]^, and we exhibit the properties of 
the 2-Toda lattice, associated with a skew-symmetric initial matrix moo(0, 0). 
The r-functions r„(t, s) then have the property 

r„(t,s) = (-l)V„(-s,-t). 

Theorem 2.1 If the initial matrix moo{0,0) is skew- symmetric, then under 
the 2-Toda flow, moo{t,s) maintains the relation 

moc{t, s) = -moo(-s, -t)^. (2.1) 

Moreover, 

h-'S,it, s) = -{Sjr\-s, -t), h-'S2it, s) = iSj)-\-s, -t), (2.2) 

/l'^^l(t, S, Z) = -%i-S, -t, Z-'), h-'^2it, S, Z) = ^t(-S, -t, Z-^), 

(2.3) 

Li{t,s) = hL^h'\-s,-t) and L2{t, s) = HLJ h'\-s, -t), (2.4) 
with h, defined by satisfying 

h{-s,-t) = -h{t,s). (2.5) 

Finally, 

r„(-s,-t) = (-l)V„(t,s). (2.6) 
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Proof: Formula ( p.l| ) is an immediate consequence of (|1.2|) and the skew- 
symmetry of moo (0,0). Formula (|2.2| ) follows from ( p.lD and the Borel de- 
composition of rriooit, s) and —m^{—s, —t)~^'- 

moo{t, s) = S^\t, s)S2(t, s), 

-moo(-s, -ty = -Sj{-s, -t)Sy^{-s, -t) 

= {Sj{-s, -t)h-\-s, -t)){-h{-s, -t)S^'^{-s, -t)), 

by the uniqueness of the Borel decomposition of rriooit, s) = —moo{—s, — t)^, 
we have 

S-\t, s) = Sj{-s, -t)h-\-s, -t) e G_ 
S^it, s) = -h{-s, -t)S^^^{-s, -t) G G+. 

Substituting {t, s) — > (— s, —t) in the second equation and comparing it 
to the first one, yields h(t,s) = —h{—s,—t), which is ( p.5|) . Substituting 
this relation into the first and second equations yields (|2.2|) , which by ( |1.5| ) 
and the definition of Li and L2, amounts to ( ^.31 ) and ( ^.41 ) . Relation ( |2.6D 
follows from (|1.8| ), ( p.l|) , footnote |^ and the multilinearity of determinant; 
or, in the semi-infinite case, from ( p.5|) , using ro(t,s) = 1: 

r„(-s, -t) Tn-i{-s,-t) ro(-t, -s) 



For a skew- symmetric initial matrix moo (0,0), relation ( |2.1| ) implies the 
skew-symmetry of moo{t, —t). Therefore the odd r-functions vanish and the 
even ones have a natural square root, the Pfaffian f2n{t)'- 

r2n+l{t, -t) = 0, T2n{t, -t) =: fl{t), (2.7) 

where the Pfaffian, together with its sign specification, is also determined by 
the formula: 



T2n{t)dXo A dXi A ■ ■ ■ A dX2n-l '■= ~] i 

^' \0<i<7<2n-l 



(2.8) 
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Theorem 2.2 For r satisfying ( \2. dj ), and hence in particular for a skew- 
symmetric initial condition moo(0,0), the 2-Toda r-function T{t,s) and the 
Pfaffians f{t) are related by 



T2n{t + [a] - [(31 -t) = f2n{t)f2n{t + N - [/?]), 
T2n+l{t + [a] - [(31 -t) = (a - (3)f2n{t - [/3])f2n+2(t +[«]), 



(2.9) 



or alternatively 



T2n{t - [(31 -t + [a]) = T2n{t - [a])f2n{t - [/5]), 
T2n(t + [a], -t - [/?]) = f2n{t + [a])f2n(t + 



r2n+l(t - [Pl -t +[a]) = {a- (3)f2n{t - [a] - [/5])f2n+2(t), 
T2n+l{t + [a], -t - [/5]) = (a - (3)T2n{t)f2n+2{t + [a] + 



(2.10) 



Proof: In formula (|1.12| ), set n = m — 1, s = — t + [/?], t' = t + [a] — [/?] and 



s' = s — [a\ — [/?] = — t — [a]; then using 



1 

27ri 



2 = 00 



f- J z = oo ^ ^ 

_^ ^ r- 1 -1 > / / r 1 1 /\ CX,Z dz 

-ReS^=p-lTm-l[t-[z \,s)Tm+l[t +[Z J,g) ^ 

(/? - «)r™_i(t - s)Tm+l{t' + S') 

(/? - «)r™_i(t - -t + [(3])Trn+l{t + [a], -t -[«]), 



/ r„(t, s - [z])rm{t\ s' + [z])e^(^-^^)^~'z"-'"-irfz 



2Txi 



2 = 

1 



(Res^=„ + Res2=/3)rm(t, s - [z\)T„,{t\ s' + [z]) 



{z — a){z — (3) 

{r^it, s - [a])Tmit', s' +[«])- r„(t, s - [(3])T^{t', s' + 



a — (3 

^-(r^(t, -t + [P] - [a])T^{t + [a] - -t) 
a — p 

- Tm{t, -t)T^{t + [a] - -t - [a] + [p])).. 
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and (|2.6|), we have 

-{(3- a)V^_i(t - 1(3], -t + [p])Tm+i{t + H, -t - [a]) 
= i-irr^t + [a] - [f3], -tf - Tmit, -t)Tmit + [a] - [(3], -t - [a] + [13]). 

Setting first m = 21 and then m = 2/ + 1, one finds respectively, since odd 
r-f unctions vanish on {s = — t} in view of 



= T2i{t + [a] - [131 -tf - T2iit, -t)T2i{t + [a] - [131 -t - [a] + [/?]), 

(2.11) 

and 

-{(3- afT2i{t - [13], -t + [l3])T2i+2{t + H, -t - [a]) 

= -T2i+i{t+[a]-[(3l-tf. (2.12) 



Taking the square root, with the consistent choice of sign,0 ( p.8| ) yields ( p^ 
and then ( [^.lUD upon setting t — > t — [a] or t — > t + [/?]. 



Corollary 2.3 Under the assumption of theorem 2.2, the wave and adjoint 



^ It suffices to check that (2.8) yields the correct sign in the second equation of (^.9| 



at = 0, < = and modulo 0{a^), i.e., 

(<9/ati)T2„+l(0,0) = f2n(0)r2„+2(0), 

for some moo (0,0) for which the right hand side does not vanish. This can be checked 
easily, e.g., for moo (0,0) made of 2 x 2 blocks i^^ 0/ diagonal. 



12 



Adler-Shiota-vanMoerbeke:Pfaff r July 23, 1999 §2, p. 13 



wave functions along the locus {s = —t} satisfy the relations 

T2n+1 



^l,2n+l(2 



s=-t 



T2n 
T2n{t) 



s=-t 



T2n 



2,2n+l\ 



T~2n-1 



^jT2n-2'T2n 



n,2n-2{< 



s=-t 



T2n-1 
T2n 



'^2,2n~liz 



-1- 



s=-t 



T2n-2 



T2r 



^2,2n-2(2: 



S = -t 



~{2n~l) 

^2n(t) 



Proof: These follow from ( |1.1(J| ), ( |2.9| ) and ( |2.1(j| ) by straightforward calcu- 
lations. ■ 



Corollary 2.4 Under the assumption of theorem \2.^ , we have 
(i) for k > 1: 

dT2n 



dtk 

dT2n+l 



s=-t 



dtk 



(a) for m > 2: 



dtkdti 

k+l=m 



:+l= 

E 

k+l=m 

^ d'^T2n 



r2n{t)^{t), 

Pk-l{-Dt)f2n ■ r2n+2{t) 

-- XI ^Pi^-^t)r2nit)){Pj{dt)T2n+2it)), 
i+j=k—l 



s=-t 



k+l=m 



dtkdti 



dtkdti 



k+l—m 



dtkdsi 



= Yl ik-l)ipki-dt)f2nmPlidt)f2n+2{t)), 
^=~t k+l=in~l 

EdT2n , s df2n / n 



s=-t 



k+l=m 
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(in) for k, I > 0: 

Pk{dt)pii-dt)T2n{t, S)\s=-t = r2nit)pk{dt)pi{-dt)f2n{t), 
Pk{dt)pi{-dt)T2n+l{t, S)\s=-t = Pl{-dt)f2n(t) ' Pk-l{dt)f2n+2{t) 

- Pl-l{-dt)f2n{t) ■ Pk{dt)T2n+2{t), 

where Pk{-) o^^e the elementary Schur functions, with P-i(-) = 0, and Dt = 
{Dtj^, (l/2)Dt2, . . . ) are Hirota's symbols. 

Proof: Relations (i) are obtained by differentiating formulas ( p.9|) in a, 
setting f3 = a and identifying the coefficients of a''~^. The first two relations 
in (ii) are obtained by differentiating formulas ( p.9|) in a and P (i.e., applying 
d'^/dadp), setting P = a and identifying the coefficients of a™-"^. The last 
relation in (ii) is obtained by differentiating the first formula in (|2.1CI|) in a 
and P, setting P = a, substituting t + [a] for t, and then identifying the 
coefficients of a™""^. Finally, expanding both identities (p.9|) in a and /?, e.g.. 



r2n{t+[a] - [P],s) = J2 c^'PW{dt)Pi{-dt)r2n{t,s) 

k,l=0 

and identifying the powers of a and P yields relations (iii). ■ 

Variants of formulas ( p.9|) and the formulas in the corollary can be ob- 
tained by using ( |2.6| ) and the following consequence of it: 

dtids-'\_, ^ ' dt-^ds'" 

where / = {ii,i2, • • • ) and J = (ji, j2, • • • ) are multiindices, \I\ = ii+i2 + - ■ ■ , 
dt^ = dfidf2 ■ ■ ■ 6tc. In particular, {d'^/dtkdsi + d"^ / dtidsk)T2n+i = 0, so we 
get the (rather trivial) counterpart of the last formula in part (ii) of the 
corollary: 

d'^'^2n+l _ Q 



dtkdsi 

k+l=m 



3 Equations satisfied by Pfaffian f-functions 

In this section, we exhibit the properties of the Pfaffian f-function introduced 
above, for the 2-Toda r-function satisfying ( |2.6| ), or the skew-symmetric ini- 
tial data moo (0,0). As in the last section, whenever we make a connection 
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with the matrix moo, we assume either moo is semi-infinite, or det m„ can be 
interpreted as in formula (*) in footnote ^ 

Theorem 3.1 The f -functions satisfy the bilinear relations 

X) 

+ / W2(t + [z\)hm{t' - [z])e^^^'''-'^^^-'z'^-'^dz = 0, (3.1) 

Jz=0 

or equivalentl'^ 



j,k>0 
j-k=-2n+2m+l 



+ Yl Pji'2y)e^~''''yk{D)f2n+2-f2n^ = 0. (3.2) 



j,k>0 
k-j=-2n+2m-l 



Proof: Formula ( ^.Ij ) follows from ( |1.12D upon replacing n by 2n and 



m 



by 2m,0 using (U) and (|2l0|) , with /3 = 0, to eliminate r2„(t - [2"^], -t), 



T"2m+i(i', -t' - [z]), T2n+iit, -t - [z]) and r2m(t' - [z], -f) and, upon dividing 
both sides by f2n{t)f2m{t)- 

d ^ i-^, i Jj, . . .), and i) = (Di, iDa, iZ^s, ■ • ■) is the corresponding Hi- 
rota symbol, i.e., P{D)f ■ g := P{d/dyi, ^d/dy2, ^d/dys, . . . )f{t + y)g{t - y)\y=Q for any 
polynomial P\ and are the elementary Schur functions: 'Yl'o' Pk{i)z^ '■— 6xp(^^ tiZ^)- 
^ One can check that all the other choices of parities of n and ra yield the same formula. 
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Substituting t — y and t + y ioi t and t', respectively, into the left hand 
side of (|3.1| ) and Taylor expanding it in we obtain 



J z=oo 

+ i e^^'y^^-'f2n+2{t - y + [z])f2m{t + y - {z])z^''-^^dz 

Jz=0 
J z=oo 
Jz=0 



z=0 



+ I ^p,(2?/)^-^eS-^»^' Y.^k{D)z''f2n^2 ■ f2 
"'2=0 j=0 fe=0 

\-fc=-2ri+2m+l 

+ 5^ P,(22/)eS-^^^>,(D)f2„+2-r2™) 

k-j=-2n+2m-l ' 



showing the equivalence of ( |3.1| ) and (|3.2| ). 



The identity (|3.1| ) gives various bilinear relations satisfied by f. We show 
that the Pfaffian f-functions satisfy identities reminiscent of the Fay and 
differential Fay identities for the KP or 2-Toda r-functions (e.g., see [0]). 
From this we deduce a sequence of Hirota bilinear equations for f, which can 
be interpreted as a recursion relation for f2„(t). 
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Theorem 3.2 The functions T2n{t) satisfy the following "Fay identity": 

t r^n (t - tlz,] - [Q]) (t - [Q]) "^-^^^-'f. 



i<i<' 

(3.3) 

i/ie "differential Fay identity" : 

{f2„(t - [M]),f2„(t - H)} 

+ iu^^ - ^^"^)(r2n(t - [U])f2nit - [v]) - f2„(t)f2„(t - [u] - [v])) 

= UV{U - V)f2n-2(t - [u] - H)f2„+2(t), (3.4) 

and Hirota bilinear equations, involving nearest neighbors: 

Pk+4{D) - ^DiDk+-^ f2n ■ f2n = Pk{D)f2n+2 " ^2n-2- (3.5) 

Here, in l \37^ ) 2n, 2m e A, I, r > such that r-l = 2n- 2m, Zi (I <i <l) 
and (i (1 < i < r) are scalar parameters near 0; in ( ^.4) 2n — 2 & A (hence 
2n, 2n + 2 G A), and u and v are scalar parameters near 0; and in ( \3.^ 
2n~2 e A, k = 0, 1, 2, ... , and {f,g} := fg - fg' = D^f ■ g %s the 
Wronskian of f and g, where ' = d/dti. 

Proof: The Fay identity ( |3.3| ) follows from the bilinear identity ( |3.1| ) by 
substitutions 

t^t- [zi] [zi] and t'^t- [Ci] [Cr]. 

Indeed, since r — I = 2n — 2m, we have 

eMY('^-t[)z')z'-'-'dz = \^'r='^^~ '";U -'-'dz 



nLi((iA)-a) 
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and 

^ nU(i-w^) 

SO the first and second terms on the left hand side of (pTTP, divided by 27rz, 
become, respectively, 

1 



- Yl ^^^^=C-' r2nlt- ^[Zj-] -[z M 

i=l ^ j=l ^ 



X^Resc=c,r2jt-5^y-[C]) 

i=l ^ j=l ^ 

i; (i - 1] 1^,1 - ic.i) f2„« (i - 101 + icj) ^ 

4=1 ^ 7 = 1 ^ ^ 7 = 1 



l<fe<r(Ci ~ Ck) ' 



and 

1 

27ri 



2=0 



Res^=^, f2n+2 f t - XI [^i] + [^] ) 

i=l ^ j=l ' 

i-1 .7-1 .7-1 
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showing 

Note that when 2m = 2n — 2, I = 1 and r = 3, denoting Zi = Q^i for 
2 < z < 4, and multiplying both sides of ( p.3| ) by Y[2<j<k<4(^j ~ ^k), we 
obtain 

(^2 - Zi){zs - Zi)f2n{t - [zi] - {Z2])f2n{t " N] " [2^4]) 

- (2:3 - 2;i)(2:2 - Zi)f2n{t - [2:1] - [zz])T2n{t - N - [2:4]) 
+ {Zi - Zi){z2 - Z'i)f2n{t - [Zi] - [Zi])f2n{t - [Z2\ - [z^]) 

+ ( n (-^^ ~ zMf2n+2{t)f2n-2{t - [Zl] - [Z2] - [z^] - [z^]) = 0. (3.6) 

The differential Fay identity ( p.4|) follows from (|3.6|) by taking a limit (set 
2:4 = 0, divide by 23 and let Zj, 0). Alternatively, we can prove ( p.4{ ) 
directly from ( p.l|) : Set t — t' = [u] — [v], 2m = 2n — 2 in ( p.l|) and in the 
clockwise integral about z = 00, set z 1/ z, yielding 



r2n{t - [Z])f2n{t' + [Z])'^ ^ ^ '^^ 







1 — u/ z z^ 

r2n+2{t + \z\)f2n-2{^ - N)^^ ^Z^dz. 



l-u/z 2 







1 — f /2; 



The first integral has a simple pole aX z = u and a double pole at z = 0, 
while the second integral has a simple pole aX z = v only, yielding, after 
substitution t' = t — [u] + [v], 

T2n{t - [u])f2n{t + {v]){u - v)^ 

+ ^ [r2n{t - [Z])f2n{t - [u] + [v] + [^])^^ 

az \ z u ^ 

= -T2n+2{t + H)r2„-2(i " - u)v^ , 

or, after carrying out d/rfzl^^o on the left hand side, 

T2n(t - [u])f2„(t + M)(U - V)^ 

V — U V 
+ T2n{t) ■ f2n{t - [u] + [v]) — - Di T2n{t) ■ hnit - [u] + [v])- 

u 

= -f2n+2(t + [v])T2n-2it - [u]){v - u)v^ . (3.7) 
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Shifting t t — [v] and multiplying both sides hj u/v yield (|3] 

Since P{—D)f ■ f = P{D)f ■ f by the definition of Hirota operator, ( |3.5|) 
is the same as (??), which, as we have pointed out, are nothing but the 



coefficients of yk+3 in ( |0.?| ), or ( |3.2| ). It also follows from ( |3.4| ), since, for any 
power series F{t,t') which satisfies F{t,t) = 0, 

f d d \ 

coefficient of yk+3 in F{t -y,t + y) = ( — - —1 F{t, t) 

= 2^F{t,t') = 2 X coefficient of u''^^ in -^F{t,t- [u] + [v]) 



Indeed, differentiating (|3^), which is equivalent to ( p.4|) , in v, setting v = u 
and using Dif ■ / = 0, 

^{Dif{t) ■ git + [v])) = -^D^DJit) ■ g{t + [v]) 

= -lj2^v^-'D,D,fit)-git+[v]l 

etc., we have 

1 11°° 

-T2n{t - [U])f2n{t + [u])^ + f2n{tf— + -^U^~^DiDjf2n " f2n{t) 
= -r2n+2{t + [u])f2n-2{t - 



which, noting f{t + [u])g{t — [u]) = ^ u^pk{D)f ■ g, is a. generating function 
for 



As in the case of KP or 2-Toda r-functions, Pfaffian r-functions satisfy 
higher degree identities: 



Theorem 3.3 

Pfaff 



/ jZj ~ Zi)f2n-2it - N - [Zj]) 
V r2n{t) 



l<i,j<2k 

T2n-2k [t - YlTl^i 

= A(^) VtT^ ^, (3.8) 



wherek > 1, 2n—2k G A, Zi are scalar parameters nearO, andA{ui, 
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Proof: This may be obtained, up to the sign, from the second identity in 
Theorem 4.2 of 



dct^ ^^'^^^' ^^'^'^^ ^^■''^^ 



l<i,j<k 



^(y)^(^) — ^ — Tins 

TN[t, S) 



by setting N hh^ 2n, k h-* 2k, Ui = Zi, taking the square roots of both sides and 
using ( |2.1CI| ). Rather than taking this route, here we prove ( p.8| ) by induction 



on k, using the bihnear Fay identity (|3.3| ). First, ( p.8|) is trivial when k = 1. 
(Note also that it gives ( |3.6| ) when k = 2.) Suppose ( p^ ) holds for — 1. 
Then we have, for every p G {2, . . . , 2k}, 



Pfaff 



f {Zj - Zi)f2n-2{t - [Zj] - [Zj]) ' 
\ r2n(t) J2<i,j<2k 



^ T2n-2k+2 (t - J22<i<2k,iytp[^i 
A{Z2, . . . , Zp, . . . , Z2k) 



r2n{t) 

Multiplying both sides by {-lY {zp - zi)T2n-2{t - [zi] - [zp]) / hnit) ■, summing 
it up for p = 2, ... , 2k, and using 

{-mZp - Z,)A{Z2, ...,Zp,..., Z2k) ^^'^ " 



Il2<i<2ki^i - ^l) n2<i<2fc(2p - Zi) 



and the identity 



2k 



Pfaff (aij)i<ij<2fc = ^(-l)^aip Pfaff (aij)2<jj<2fc V(ajj) skew symmetric 



p=2 ''^^P 



21 



Adler-Shiota-vanMoerbeke:Pfaff r July 23, 1999 §4, p. 22 



which follows from definition (p.8|) of the Pfaffian, we have 
Pfaff I 



./ {Zj - Zi)T2n-2it - N - \Zj\) 

V r2n{t) 



l<i,j<2k 



1 V c 



• T2n^2{t - [Zl] - [Zp])f2n^2k+2 {t- ^ [Zi] j 



2<i<2k 



using the bilinear Fay identity (|3.3|) with r = 2A; — 1, / = 1, C,i := Zi+i 
{1 < i < 2k — 1) and (2n, 2m) replaced by (2n — 2, 2n — 2A;) this becomes 

A(z) 1 



n2<i<2fc(2;i - Zl) f2n{tf 

(2k N ^ 2A; X 

W{zi - Zi) j T2n{t)f2n-2k { t - ^[Zi] j 
i=2 ' ^ i=\ ' 

T2n-2k[t - Z]^=l[^i. 



r2n{t) 

completing the proof of (|3.8| ) by induction. 



4 Vertex operators for Pfaffian f-functions 

In terms of the operators 

X{t,X) ■= e^i *fe^'e~^i 

acting on functions f{t) of t = (ti,t2, ■ ■ ■) & define the following four 
operators^ acting on column vectors g = {gn(t))n£A, 

Xi(/i):=X(t,^)x(/i), X^A) :=-x*(A)X*(t,A), 

X2ifi) := -X(s,/i)x*(/i)A, X;(A) := A^xWX*{s, A), 

^ Here X(s, A) has Si in place of ti, as well as d/dsi in place of d/dti, in the definition 
of X(t, A), etc.; x(a*) := (A*")neA, and x*(a*) = x(a*"^)- 
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and their compositions^ 



%(/i,A) :=X;(A)X,(/i) 



They form a set of four vertex operators associated with the 2-Toda lattice. 
Among those, X12 is important in the semi-infinite case, related to the study 
of orthogonal polynomials. In we showed that 



for any n, m & A, n > m. Note on the right hand side the limit exists as s — >■ 
—t if n and m are even, so in particular, taking n = m + 1, we see the poles 

along S = -t cancel out in ^l,2m(Ai)^2,2m(-^~^) + ^l,2»n+l(yU)*2,2m+l(-^"^)- 

We shall come back to this point after proving the following theorem and its 
corollary. 

Suppose r satisfies (|2.6|), and let f be the vector of corresponding Pfaffian 
f-functions. Let Xi, XI and Xn act on f as if they are acting on the extended 
vector {fn)neA, where f„ = if n is odd, so that (resp. X*W) ^-cts on 

by multiplication of /z^" (resp. A~^"). Then we have^ 

When i = j, X* interacts with Xi nontrivially, yielding the factor exp(^(/i/A)'^/A;) = 
1/(1 — /i/A) if we bring the multiphcation operators to the left and the shift operators in 
t or s to the right. So if we denote by : : the usual normal ordering of operators in s 
(but not in the discrete index n), we have 



where u — t, 4){ij., A) — x(/i/A) if i = 1; u = s, 0(/i, A) — AJ x{^l lA-^ \i i — 2] and 




(4.1) 



X,,(Ai, A) = 1/(1 - /VA)) :X*(A)X,(^): 



^l|X)cj>{^i,X)X{u,^Ji,\), 




-k 



-k 



a 



11 



The product Xi(A)Xi(/x) in (4.4) is computed the same way as in footnote |T^. 
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Theorem 4.1 

(X„(,.A)r).|.-_. = ^-W^"<''-^'^-« <'' = 2"> (4.2) 

\-A(Xi(f,)f2„(*))X;(A)?,„+2 (/V = 2n + 1) 

(X,,(,.A)r)„|.-_, = |-^-«^"(^'''>*-'*' <^ = 2"' (4.3) 

' \-fi(Xi(A)f,„)(X;(fi)f2„+2) (iV = 2n + l) 

(X.,(,.A)r)„|,_.. = |-^*-'*'^'(^)^'<'')*--'*' <^ = 2") (4.4) 
' ' ' ' \(Xi(,.)f2,.)(X;(A)f2.) (iV = 2n+l) ^ ' 

Corollary 4.2 For k = 1, 2, the following holds: 

jt\s)r2n{t,s)U=^t = {-lff2n{t)jl'\t)f^n{t), 

and so 

{J\'\t) + (-l)'=jf(,))r2„(t, = 2f2.(t)4'^(t)f2n(t). 



Proof: The theorem follows from formulas (|2.9|) and ( |2.10|) by straightfor- 
ward calculations: 



(Xii(/i, A)r)Ar| 



/ A — U 



-A/ A — 

for = 2n: 



A 

-A/ A — /i 

7'2n(t)Xii(/i, A)f2„(t) 



s AT \ 



for TV = 2n + 1: 

A(Xi(/i)f2n)(Xt(A)W2); 
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\s=-t 



{X22{IjL,X)t)n\ 

- - i-Y'^ y^e-S*^(-^-^\^(t, -t - + [A-]) 

\H/ A — jl 

\I1J A — jl 

for N = 2n: 

- - j^e-^'^^^-'%^{t)f2r.{t + - [A-]) 

= -T2n(^)Xii(A,//)f2n(i), 



for N = 2n + l: 

2" _ \-l 



\/^/ A — // 

■\ 2n 



^2n+l 



= -/x(Xi(i, A)f2n)(Xi(i,//)f2n+2); 
(Xi2(/X, A)r)7v|^=_t 

= (ATx(M)^*(5,A)X(t,/i)r)Ar|,=_i 

= (/xA)^-ieS(*'''*-(-*')^^)r^_,(t - [/x-i], -t + [A'^]) 



for N = 2n: 



= (/xA)2"-ieS*>(^'+^')(A-i - /x-i)f2„_2(t - [A-^] - [f^-'])f2n{t) 

= -A(Xi(A)Xi(/x)f2n-2(t))f2n(i), 

for AT = 2n + 1: 

= (Xi(//)f2„)(Xi(A)f2n). 
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The corollary is shown by expanding Xn in A and fi — X. Recall that 



Xii(/i, A) 



(9"a-(..a)) 



X — fi 

X f^{f.-xy 



X — jjL 



fc=0 ■ / = -oo ^ "G^ 

where A) is the vertex operator in the KP theory [P], and 

k 



(fc-j) 



i=0 

with Vr/'^'' the coefficients of similar expansion of A). 
Expanding Xn in ([4 .21 ) as above leads to 



In particular, since jf'^'* {k < 2) and (fc < 2) are linear combinations 

of each other : 

W^? = = ko, 

W^iJ = >/f + (2n - ^ - 1) J« + nin - 1) , 
one sees for k = 1, 2 that 



Consider the following vertex operator[^ 



As before, X treats f as an extended vector (f„)„gyi, where -r„ = for n odd. So x(z) 
appearing in X(z) acts on t„ as multiplication by z", and acts on f as (A^f )„ = t„-i. 
In practice, we always have even number of X's acting on f, so there is always an even 
power of A^, and Tn for odd n will never appear in our formulas. 



26 



Adler-Shiota-vanMoerbeke:Pfaff r July 23, 1999 §4, p. 27 



and define the kernel 



2ra 



It is easy to see that (X(y)X(z)r)2n = 2/Xi(?/)Xi(2;)r2„_2, so by ( ^73| 



(Xi2(^,l/)r)2n 



T2n 



(Xi2(y,^)r)2„ 



s=-t 



T2r: 



and by (g^) 



2m<j<2n 



Knifl, A) - Kmifi, A). 



s=-i 



Here each term \l/ij(/i)\l/2 j(A~^) on the left hand side blows up along s = —t, 
but the poles from two successive terms (for j = 2k and j = 2k + 1) cancel, 
as we saw earlier. 



For n E A, let 



' 2m r2m{t - [)^-']) 



qn{t,X) := < 



^/r2m{t)r2m+2{t) 

^^(d/dh + x)f2Ut-[x-']) , 



if n = 2m, 
if n = 2m + 1. 



(4.5) 



a/ 'T~2rn{t)'T2m+2{t) 

In the semi-infinite case, the g^'s form a system of skew-orthogonal polyno- 
mials [0. 

Theorem 4.3 The following holds: 

Vm{Kr,{z,,z,)\<,,,<2k=i\ \{MZi)A , (4.6) 

\ «=1 , / 2n 



Kn+l{fi, A) - Kn{^, A) 

so m i/ie semi-infinite case 



ordered 



gEr*.(M'+A0(^2n(t,A)g2n+l(t,/^) -g2n(t,/i)g2n+l(t,A)), (4.7) 



AT-l 



KNifl, A) = e^^*'('^'+^') {q2n{t, X)q2n+l{t, fi) - g2n(t, fi)q2r^+l{t, A)) 



(4.8) 
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Proof: Using (|3.8| ) and 

/X(^)X(A)f 



2n 



(A-i-/i-^)(/iA)2"-ie^ 



r f,(A^'+A') ^2n-2(t - [/X ^] - [A ^]) 



(4.9) 



the left hand side of 



becomes 



/ \2n-l T'^'' T°° t z^ \ /- _ix^2n-2fe(i J2j=l[^j ]) 



T-2n(t) 



This equals the right hand side of (WM), because 



2n 



2n-2 ^2n-2 2n-2fc J] ti{z{+--+zl ) 



n n(i 

•l<jr<2fc l<i<j 



T2n-2fc ^ 



1 / 

(zi ■ ■ ■ Z2kr-'A{z-')e^ *'(^^+-+^-)f2„-2. (t - f;[^r']) • 



To prove (|4.7|), we have 



(/i - A)(/iA) 



T"2n+2(i) T2n(i) 
2„T2„(t)f2„(t - [f^-'] - [A-1]) - (//A)-2f2„-2(t - [yW-^] - [A-l])f2n+2(t) 



(/iA) 



T'2n+2(i)T-2n(^) 

2„{r2„(t - f2„(t - [A-1])} + (/X - A)f2n(t - [/i-'])r2n(t - [A^^]) 



T2n+2(i^)T2n(^) 



using (13:^), 



,2n r2„(t-M) /i'"(5/9ti+/x)f2n(t-[/i-^]) 



^/T2n{t)f2n+2{t) \/ T2n{t)T2n+2{t) 

(g2n(^, A)g2n+l(i,/i) - g2n(i, /i)g2n+l (i, A)) 



(A ^ /i) 
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in terms of the skew-orthogonal polynomials ([4.5|). Multiplying this with an 



exponential and noting ([4 .91) , one obtains (|4.7| ). Summing up this telescoping 
sum yields 



X(/i)X(A)f 



' 2N 



T2N[t) 



N-1 



J2 e>^^'^^^'+^'\q2n{t, X)q2n+i{t, fi) - q2n{t, fi)q2n+i{t, A)). 







5 The exponential of the vertex operator main- 
tains f-functions 

The purpose of this section is to show the following theorem: 
Theorem 5.1 For a Pfaffian f -function, 

f + aX(A)X(/i)f (5.1) 
is again a Pfaffian f -function. 

Remember that X(A)X(/i) acts on f2„(t), as follows: 

X{X)X{fi)f2n{t) = (l - ^) A2"-V"-'eS*^(^'+^')f2„_2(t- [A-^] - [fi-']). 

(5.2) 

It is convenient to relabel —>■ Tn 

X(A)X(/.)f„(t) = (l - 9 A^"-V"~^eS?°*'(^^+''^)f„_i(t - [A-i] - [/.-I]) 



^(A-/i)Y(A,/i). (5.3) 

29 



Adler-Shiota-vanMoerbeke:Pfaff r July 23, 1999 §5, p. 30 
With this relabehng, the bihnear identity takes on the form 



+ i Tn+iit + \z])T^{t' - [z])e^^^''^-'^>"z^^-'"'dz = 0. (5.4) 



2=0 



Lemma 5.2 We have: 



1 .^r, iv^^r.-ii. 



/i — Ay Ay A — /x \ /i^ 

Proof: See for instance 
Lemma 5.3 

+ / Yr„+i(t + [z])TUt' - [z])e^^^<~'^>-\'--'"^dz 



z=0 



Proof: Upon performing the following operations 
n I— > n — 1 

multiphcation by (A/i)^"~^e^i° 



(5.5) 
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the bilinear identity (5.3) yields 

J z=oo 

^1 _ ^1 _ ^ eEr((*-«i)^*+*i(M'+A'))^2n-2m-2^^ 



2n 



+ i r„(t + [z] - [/i-i] - [A-^])r„(t' + [^])(A/i) 
Subtracting this expression (which is = 0), the left hand side of (5.5) equals 

/ Tn-l{t - [Z-'] - [A-1] - [/i-^])wi(t' + [Z-']) 
J z=oo 

gEr((t.-t^)^'+(t.-^)(A'+M0)(;^^)2(n-l)^2n-2m-2^^ 

+ / Tn{t + [Z] - [\-'] - [^-^])Tm{t' - [z]) 
Jz=0 

gEr((*^*.)^-^+(*.+4)(A'+M0)(;^^)2„^2n-2.„^^ 

= / T^-i{t - [z-^] - [X-'] - [i^-'])T^+,it' + [^-l])eI:^((*.-*^y-^*.(A^+MO) 

J z=oo 

+ / r„(( + Iz] - [A-'l - [M-'l)r„(i' - W)eS"(«-'*-+"(*'+'''»(A„) 
'd - A.)-'(l - - ^ (l - ±) " (l - ± 



2n 



1 



n — X 



- A^>^-r„(t - [A-^])r^(t' - [/.-^])eS^°(*^^^+*>^) | , 
ending the proof of the lemma. 
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Proof of theorem 5.1: It suffices to prove 

J 2 = 00 

Jz=0 

The coefficient of and 6^ vanishes, on view of the fact that r„ and Yr^ are 
Pfaffian r-functions. So it suffices to show the vanishing of the a6-term. 

afe-coefficient 

= (f (YT„(t - [z-'])rr^+,{t' + [Z-']) + Tn{t - [z-'])YTrn+,{t' + [z"^])) 
J z=oo 

eEr(*i-*'i)^*^2n-2m-2^^ 

+ / (YT„+i(t + [^])r^(t' - [^]) + r„+i(t + [;.])Yr^(t' - [^])) 

Jz=0 

^2n-2m^^_ 

The first terms in each of the integrals can be evaluated by means of lemma. 
The sum of the two terms equals 



1 



jji — \ 

-A^>^™r„(t - [\-'])Tm{t' - [/x-^])eS?°(*'^'+*>^) j . (5.6) 

Performing the exchange 

n < — > m, t < — > t', z < — > z~^ 

gives an expression for the sum of the second terms in the integrals; the 
sum of expression (5.6) and the same expression with the exchange above is 
obviously zero. ■ 

6 Examples 
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6.1 Symmetric matrix integrals 

Consider the matrix m„(t, s) of {t, s)-dependent moments, 

fiUt,s):= [ [ x''y'e^^^''^'-'^y'^Fix,y)dxdy, t,seC'^, (6.1) 



with regard to a skew- symmetric weight F{x,y), satisfying 

F{x,y) = -F{y,x). 

Then 

Tn{t, s) := det m„(t, s) (6.2) 
= /■•■/ n(eS-i(*»^'^-^>^^-)F(xfe,yfc)) A„(x)A„(y)rfi4 (6.3) 



The proof that the determinant (|6.2|) of the moment matrix equals the 
multiple integral ( |6.3| ) is based on an identity involving Vandermonde deter- 
minants and can be found in 0. Then 

f2„(t) := Pfaff m2„(t, -t) = A/r2„(t, -t) 

satisfies equations (0.6) up to (0.8). 

Moreover, the moments fiij in (|6.1| ) satisfy the equations 



'W, — ~ ^^^+k,j ana -- — — — /ij 
otk dsk 



and so m := moo satisfies (|1.1|). 

The skew-symmetry of F above implies the skewness of moo(0, 0); so, by 
theorem 2.1, we have 

/iij(t, S) = -IJ.ji{-S, -t). 

Therefore also, moo, Si, 5*2, ^^i, ^^2, h and r have the properties mentioned 
in theorem 2.1 and the vector (x„)„>o is a 2-Toda r-vector, i.e. the corre- 
sponding wave vectors and the matrices Li satisfy the formulae of theorem 
2.1. 

For skew-symmetric weights F{x,y) of the special form 

F{x,y) := e^^^^'^'^^'^'^ lE{x)lE{y) sign(a; — y), for an interval _E C M, 
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and for a union of intervals C M, the expression T2n = T2n(i, — t)^^^ equals 
the integral over symmetric matrices, given in 

for the Haar measure dX on symmetric matrices and 

S2n{E) := {2n x 2n symmetric matrices X with spectrum C E}. 

In 01, we worked out the Virasoro constraints satisfied by (??), which 
then leads to inductive expressions for those integrals, involving Painleve- 
like expressions. 



6.2 Quasiperiodic solutions 

In this subsection, we shall combine the construction of quasi-periodic so- 
lutions of 2-Toda lattice |l^, and the theory of Prym varieties [|l5l to 



obtain quasiperiodic solutions of the Pfaff lattice. While we put stress on 
the semi-infinite case in the present paper, this gives a non-trivial example 
in the bi-infinite case. 

A 2-Toda quasiperiodic solution is given by some deformation of a line 
bundle £ on a complex curve (Riemann surface) C, with the time variables 
playing the role of deformation parameters, so the orbit under the 2-Toda 
flows is parametrized by the Jacobian of C. If C is equipped with an invo- 
lution l: C ^ C, and if £ satisfles a suitable antisymmetry condition with 
respect to l, then the 2-Toda flows can be restricted to preserve the antisym- 
metry of C, giving a solution of Pfaff lattice. The Prym variety P of (C, l) 
naturally appears as the restricted parameter space. The vanishing of every 
other Tn{t, —t) (see ( |0.4| ) or (p^.6D) indicates that the space of £'s which satisfy 
the antisymmetry condition must consist of two connected components, P 
and P~ . This means the involution t has no flxed points. So, in general a 
quasiperiodic solution of the Pfaff lattice does not satisfy the BKP equation 
and vice versa, since the orbit of a quasiperiodic solution of the BKP equa- 
tion is isomorphic to the Prym variety of a curve with involution having at 
least two fixed points. 



Preliminary on the geometry of curves 

A line bundle on a complex curve C is defined by a divisor D = f^iPi 
rrii E Z, Pi E C, i.e., a set of points Pi on C with (positive or negative) 
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multiplicities m^, as £ = 0{D). Its local sections (on an open set U (Z C, 
say) arc mcromorphic functions on U which have poles of order at most 
rrii (zeros of order at least —rrii) at Pi. The number d := 'Y^nii is called 
the degree of C For C = 0{D) and m, n E Z, p, q E C, we denote 
jC{mp + nq) — 0{D + mp + nq) etc. A deformation of C can be described 
as a deformation of D, hke Dt s = Yl/T^iPiit, s), but in the 2-Toda theory it 
is more convenient to describe it by requiring its local sections to have some 
exponential behaviors at prescribed points, as we shall see later. 

Two line bundles 0{Di) and 0{D2) are isomorphic if the divisors Di 
and D2 are "linearly equivalent," i.e., if they differ by the divisor of a global 
meromorphic function on C. Jacobian J of C is the space (Lie group) of 
isomorphism classes of degree line bundles on C . It becomes a principally 
polarized abelian variety of dimension g := genus of C, i.e., J is a complex 
torus C^/r, D r ~ Z^^, for which there is a divisor (codimension 1 
subvariety) © C J, such that some positive integer multiple of © defines 
an embedding of J into a complex projective space, and © is "rigid" in the 
sense that it has no deformation in J except parallel translations. A complex 
torus O/r is a principally polarized abelian variety if and only if, after some 
change of coordinates by GL{g, C), the lattice F becomes + QZ^ for some 
complex symmetric g x g matrix Q with positive definite imaginary part. 
On a principally polarized abelian variety C^/(Z^ + QZ^), there is a special 
quasiperiodic function (i.e., holomorphic function on that satisfies some 
quasiperiodicity condition with respect to Z^ + ilZ^) called Riemann's theta 
function t?, defined by 



The theta divisor © becomes the zero divisor of 

If C has a (holomorphic) involution l: C ^ C (i.e., = id), J gets an 
involution l* induced by l. The Jacobian J' of the quotient curve C = C/l, 
and the Prym variety P of the pair (C, c) (or (C, C')) appear in J roughly as 
the ±1 eigenspaces of l: J' — J'/ (some subgroup of order 2) C J and P C J 
are subabelian varieties of J, such that ijj, = +1, l\p — —1, and J ~ (J' x 
P)/(finite subgroup). When l has at most two fixed points, the restriction 
of on P gives twice some principal polarization on P (the restriction ^9|p 
becomes the square of the Riemann theta function on P defined by this 
polarization) . 
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Quasiperiodic solutions of 2-Toda lattice 

Let C be a nonsingular complete curve on C (compact Riemann surface) of 
genus g, let £ be a line bundle of degree — 1 on C, let p, g G C be distinct 
points. Let us choose local coordinates at p and z at g, and trivializations 
of C{j)) at p and g, 

dp: Cp{p) ~ Op and a^: CqC:i Og. 

For t, s G C°°, let Ct^s be the line bundle whose (local holomorphic) sec- 
tions are (local holomorphic) sections of C away from p and q, and at p 
(resp. q) have singularities of the form ■ (holomorphic) (resp. '^'^ ' • 

(holomorphic)). For "generic"Q {n,t,s) G Z x C°° x C°°, the wave functions 
^i,n5 ^2,n are obtained from a (unique) section ipn(t,s) of 

£i_s((n + l)p-ng), 

which has the form z"'e^*'^*(l + 0(2;"-'^)) at p via CTp, i.e., 

v[/i,„(t, ^) := ap(^„(t, s)) = z-e^'^^\l + Oiz-')), 
^2,n(t, s, z) ■= ag((^„(t, s)) = z^e^ '^''\K{t, s) + 0{z)). 
The adjoint wave functions 

^tn = -2""e-^*^^'(l + 0(;2-i)), 

n,n = ^-''e-^'''~\hnit,s)-' + Oiz)) 

are defined similarly, by using 

{Ct^sTi-np +{n + l)q) = {C*)^t,-s{-np + (n + l)g), 

in place of Ct,s{{n + l)p — nq), where we denote 

£* := Hom{C, u) = ® u, 

with uj being the dualizing sheaf (the canonical bundle, i.e., the line bundle 
of holomorphic 1-forms), and, in place of o"p and a^, trivializations 

a;:C;c:^O p and a;-. C;{q) c:^ Og, 

Here generic means that r{Ct,s{np — nq)) ~ {0} holds. For a degree g — 1 hne 
bundle £, this condition holds for almost all {n,t,s) e Z x C°° x C°°, and implies that 
dimr(£t^s((n + l)p - nq)) = 1. 



36 



Adler-Shiota-vanMoerbeke:Pfaff r July 23, 1999 §6, p. 37 

for which the maps 

Cp{p)(^Cl 3 ^ (7p{(f))a*{tlj)dz/z e uj{p)p, 

® C*g{q) 3 (0, ij) ^ ag{(p)(r*{ij)dz/z G uj{q)g 

extend to the canonical map 

C{p)0C*{q) ^ uj{p + q). 

Hence for general (n, t, s), (m, t', s') G Z x C°° x C°°, 

^i,„(t, s, z)^l^{t', s', z)dz/z, 2 = 1, 2 

become expansions at p and q, respectively, of a holomorphic 1-form on 
C \ {p,q}, so by the residue calculus the pair satisfies the bilinear 

identities (|1.7|). 

Quasiperiodic solutions of Pfaff lattice 

In the above construction, suppose C has an involution l: C ^ C with no 
fixed point. In this case g is odd, g = 2g' — 1, with g' being the genus of the 
quotient curve C = C/t. Suppose q = l{p), and C satisfies 

L*{C)c^C*, so that C0i*C-iu. (6.6) 

Choose the local coordinates z"^^ and the trivializations cxp, ag, a*, a* at p 
and q = t{p), such that z ■ l*z = 1 and ag = l* o a* o l* hold. (We then 
have a* = — i* o dp o 6*, with the minus sign due to the fact that dz/z, which 
appear in ( |6.5| ), satisfy L*{dz/z) = —dz/z.) Then the wave and adjoint wave 
functions constructed above satisfy (|2.3| ), so they lead to a quasiperiodic 
solution of the Pfaff lattice when s = —t (and skipping every other n). 

The orbit of the 2-Toda flows is parametrized by the Jacobian J of C, 
and the r-functions are written in terms of Riemann's theta function of J . 
The orbit of the Pfaff flows will become the Prym variety P of (C, t), with f 
given by the Prym theta function. To be more precise, let Jg-i be the moduli 
space of the isomorphism classes of line bundles of degree g — 1 on C . This 
is a principal homogeneous spaceQ over J, on which the theta divisor 

©:={/ : G J,_i I r(£) ^ (0)} 

Hence Jg-i is (non-canonically) isomorphic to J. We choose this isomorphism in such 
a way that Q C Jg-i is identified with the zero locus of Riemann's theta function for J. 
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is canonically defined. The set of £ G Jg-i satisfying ( |6.6|) becomes the 
disjoint union Pg-i U Pg^i, where 

Pg-i := G Jg-i I C satisfies and dimr(£) is even}, 
Pg-i '■= G Jg-i I C satisfies and dimr(£) is odd} 

are principal homogeneous spaces over the Prym P. We have 

Pgli C e and Pg_i ■ = 2H, 

for some divisor S C Pg~i which gives a principal polarization on Pg-i- 
Since is the zero locus of Riemann's theta function d of the Jacobian J, 
this means d vanishes identically on Pg^i-, and the restriction ^\pg_^ becomes 
the square of Riemann's theta function dp oi (P, S), which is called the Prym 
theta function. 

For a 2-Toda quasiperiodic solution, the discrete time fiow (shift of n by 
1) is given by the shift C ^ C{jp — q) . In the present case, since q = tip), 
this fiow preserves condition (|6.6|) . Moreover, we have 



yp ec,yce Jg^, . I ^ ^ ^ _ ^^^^^ ^ (6.7) 

so that C{np — m{p)ys alternate between Pg-i and Pg_i, and every other r 
function vanishes identically when s = —t. Shifting the discrete index n by 



1 if necessary, we may assume that Tn{t,s) satisfies (|0.4|) or {^M 



Explicit formulas 

Explicit formulas for \E', \E'* and r can be given in terms of Riemann's theta 
function for J, and hence explicit formulas for f can be given in terms of the 
Prym theta function for P. 

Taking a basis Ai, Bi {i = 1, . . . ,g) of Hi{C, Z) such that Ai ■ Bj = 6ij 
and Ai ■ Aj = Bi ■ Bj = 0, let Ui {i = 1, . . . ,g) be a basis of the space of 
holomorphic 1-forms such that 



Then 
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gives a complex symmetric matrix Q with positive definite imaginary part, 
and J = C^/ (Z^ + becomes the Jacobian of C. Choosing a point p & C, 
the map 

a: C 3 X ^ (^j ui, . . . , J tu^ E J 

is well-defined and gives an embedding of C into J. Composing a with a 
translate of Riemann's theta function: 

q{x) ■.= 'd{a{x) + a), a G (6.8) 

one obtains a multi-valued function on C which is single-valued around the 
y4-cycles. 

Next, let (n \ n = 1, 2, . . . , be the differentials of the second kind 
(meromorphic 1-forms with no residues) with poles only at p of the form 
d(z" + 0(1)) and no A-periods (/^^ C"^^ = 0), and let ci^\ n = 1, 2, . . . , be 
defined similarly, with p replaced by q and z by (recall that (resp. z) 
is the local coordinate at p (resp. q)). Let Co be the differential of the third 
kind (meromorphic 1-form with simple poles) with no A-periods and poles 
only at p and q of the form dz/z + 0{l). Then, given (n, t, s) G Zx C°° x C°°, 
the multi-valued holomorphic function 



C 3 X e(x] 



/ UC0 + 5^U^^) + 5^^nC^'^M (6.9) 
^ n=l n=l ^ ^ 



has singularities at p and q of the form z'^e^ and z"'e^ respectively, 
and is single- valued around A-cycles. The product of the form e{x)q{x) / q{p), 
where e{x) and q{x) are as in (|6.8|) and (|6.9|), with 



a = a(n, t, s) = nQ;(g) + + ^SjVi + Oo, Voq G C^, (6.10) 



i=l i=l 



and t/j = —{d/d{z ^)ya{p)/{i — l)\,Vi = —((i/(i2;)*a(g)/(2—l)!, gives function 
on {n,t, s), and hence the wave functions \E', with the desired properties: 

• ipnit, s] x) is single- valued around the A-cycles, and when x goes around 
Bi, it is multiplied by a factor independent of {n,t, s), 

• v^„(t, s; x) ~ z"'e^**^*(l + 0{z^^)) at x ~ p, and 
(Pnit, s', x) ~ z'^e^ ' ihn{t^ s) + 0{z)) at a; ~ q. 



39 



Adler-Shiota-van Moerbeke : Pf af f r 



July 23, 1999 §6, p. 40 



The adjoint wave functions ^* are obtained similarly, from e{x)~^'3{a{x) — 
a)/'d{—a) with the same a as above. 

The 2-Toda r-function can be computed from those formulas as 

Tn{t, s) = exp((5(n, t, s))-&{a{n, t, s)) (6-11) 

for some quadratic form Q{n, t, s), i.e., 

oo oo oo 

Q{n, t,s) = ^ QijUtj + Q'ijSiSj + n{qiti + g-s^), 

i,j=l i,j=i i=l 

with Qi j = Qj i appearing in the Laurent expansion of the integral of (j^^ or 
) as ' ' 

r = - 2 g,,,.^-Vj for xc,p, 

Q'i. j = Q'j,i appearing similarly in the Laurent expansion of the integral of 

C^^ or as 

and qi and q[ appearing similarly in the expansions 



Co = log 2; - ^ qjZ"^ I j for 

.7=1 



a; ~ j? 



and ^ 

/ Co = logz — QjZ'' /j for a; ~ g. 

Suppose C has an involution l with no fixed points, so that g = 2g' — 1 
with g' being the genus of the quotient curve C — C/i. Suppose q — i{p). 
Take the cycles Ai, {i = l,...,g) in such a way that L{Ai) ~ 
L{Bi) ~ Then 6*(ci;j) = u;g+i_j, and ^2 satisfies Qij = Qg+i^i^g+i^j. 

The map l: O 9 . . . ,Zg) i— >• . . . , 2:1) e O maps the lattice F := 

+ QZ^ onto itself, and the embeddings 

J' = J'/{Z/2Z) C J and P C J 
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are given by the images under tt^ : /T of the ±l-eigenspaces of I: 

Setting 



R' '■— ~^ ^i,g+l~j)l<i<g,l<j<g' R" ■— {^ij ~ ^i,g+l- 



-j)i<i<g,i<j<g'-i^ 



so that := R'C^' and := R"C^' ^ are the ±l-eigenspaces of l, and for 
any z G z' := {l/2){R'Yz and z" := (l/2)(_R")*z give the decomposition 



z = R'z' + R"z" G CI © C^, we have 



z' ^ R'z' 

and 

P = a'-^/{Z<^'-^ + Q"Z9'-i) ^ 7ri(i?"Cf'-^) C CV(Zf + Q'Z^ 



where e = diag(l, 1, . . . , 1, 1/2), 



i<i,j<g' 



In ( |6.10|) , suppose oq = R"a'Q G -R"C^' ^. Since, by definition, = a(g) 



a{p) G 7ri(Cl) and = V^, we then have a{n,t, —t) = R"a"{n,t), where 

1 1 °° 

a"(n,t) = -{R"ya{n,t,~t) = {R")' (^-na{q) +J2t^U?) + a'^. 



1=1 



Hence by using ( |6.11D and ( p.l2| ), and noting that Q- = Qij and q'^ = —qi 
we have 

fit) = exp{Q{n,t))'dp{a"{n,t)), 

where 

i,j = l 1=1 

and 

■(9p(z) = ^ exp{2']Ti'm^z + 7iim^fl"m), for 2; G C^'"^ 
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